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The  kinetics  of  the  reaction  of  with  O2  has  been  studied  from  120  to  1400  K  using  both  a 
selected  ion  flow  tube  and  high-temperature  flowing  afterglow.  The  rate  constant  decreases  from  120 
K  to  ~1200  K  and  then  increases  slightly  up  to  the  maximum  temperature  studied,  1400  K.  The  rate 
constant  compares  well  to  most  of  the  previous  measurements  in  the  overlapping  temperature  range. 
Comparing  the  results  to  drift  tube  data  shows  that  there  is  not  a  large  difference  between  increasing 
the  translational  energy  available  for  reaction  and  increasing  the  internal  energy  of  the  reactants  over 
much  of  the  range,  i.e.,  all  types  of  energies  dinve  the  reactivity  equally.  The  reaction  produces  both 
NO'*‘  and  NO^  ,  the  latter  of  which  is  shown  to  be  the  higher  energy  NOO'*'  linear  isomer.  The  ratio 
of  NOO'*'  to  NO'^  decreases  from  a  value  of  over  2  at  120  K  to  less  than  0.01  at  1400  K  because 
of  dissociation  of  NOO"^  at  the  higher  temperatures.  This  ratio  decreases  exponentially  with 
increasing  temperature.  High-level  theoretical  calculations  have  also  been  performed  to  compliment 
the  data.  Calculations  using  multi-reference  configuration  interaction  theory  at  the  MRCISD(Q)/ 
cc-pVTZ  level  of  theory  show  that  singlet  NOO'*'  is  linear  and  is  4.5  eV  higher  in  energy  than 
ONO'*’.  A  barrier  of  0.9  eV  prevents  dissociation  into  NO'*'  and  0(’Z));  however,  a  crossing  to  a 
triplet  surface  connects  to  NO'*'  and  OCP)  products.  A  singlet  and  a  triplet  potential  energy  surface 
leading  to  products  have  been  determined  using  coupled  cluster  theory  at  the  CCSD(T)/ 
aug-cc-pVQZ  level  on  structures  optimized  at  the  Becke3-Lee,  Yang,  and  Parr  {B3LYP)/ 
aug-cc-pVTZ  level  of  theory.  The  experimental  results  and  reaction  mechanism  are  evaluated  using 
these  surfaces.  ©  2004  American  Institute  of  Physics.  [DOI:  10.1063/1.1807376] 


I.  INTRODUCTION 

Air  plasmas  at  elevated  pressures  are  currently  the  sub¬ 
ject  of  intensive  study  for  a  variety  of  applications.  Typical 
nonequilibrium  plasmas  have  heavy  particle  temperatures 
below  1000  K  and  average  electron  energies  on  the  order  of 
a  few  eV.  These  conditions  are  favorable  for  associative  ion¬ 
ization  processes  and  the  presence  of  polyatomic  nitrogen 
ions  such  as  and  N4  is  expected.  However,  temperature 
dependent  data  even  for  the  relatively  simple  cases  involving 
reactions  of  these  ions  with  major  species  are  insufficient  or 
nonexistent. 

In  light  of  these  issues,  the  kinetics  of  the  reaction  of 
with  O2  has  been  studied  in  more  detail.  The  reaction  has 
two  observed  product  channels: 


^N0+-F(N20)  +  466  kJmoF* 

(la) 

^N02+N2+124  kJmoF’, 

(lb) 
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where  the  N2O  product  may  be  dissociated,  i.e.,  N2  +  O,  re¬ 
ducing  the  exothermicity  of  the  reaction  to  299  kJmoF*,  or 
production  of  N+NO,  which  would  make  the  reaction  endo¬ 
thermic  by  14  kJmoF'.  The  heat  of  formation  of  has 
been  obtained  by  combining  the  average  of  the  experimental 
heats  of  fonnation  of  N3  from  the  work  of  Continetti  et  al.^ 
and  Martin  et  al}  with  the  ionization  potential  from  the  work 
of  Dyke  et  al.^  which  are  in  good  agreement  with  the  recent 
theoretical  values  of  Dixon  et  al.^  For  the  NO^  thermochem¬ 
istry,  we  use  the  present  ab  initio  calculations  for  the  NOO'*' 
isomer  (see  below).  All  of  the  other  values  are  from  the  NIST 
Webbook.^  Note  also  that  the  NIST  Webbook  recommended 
value  for  the  heat  of  formation  of  N3  is  52  kJmoF'  lower 
than  the  values  used  here. 

In  a  flowing  afterglow,  Dunkin  et  al.^  measured  the  first 
experimental  value  for  the  total  rate  constant  for  reaction  (1) 
at  200  K  and  obtained  a  value  of  k=  1  X  10“*®cm^  s“'. 
Smith  et  al.^  have  determined  a  rate  constant  of  5.1 
X  10~^'  cm^  s'"’  for  the  reaction  at  300  K  in  a  selected  ion 
flow  tube  (SIFT),  which  is  an  order  of  magnitude  less  than 
the  Langevin  collision  rate  constant  of  7.0X  10“'°cm^  s“'. 
In  addition,  Lindinger*  has  measured  the  rate  constants  for 
the  N^  reaction  with  O2  in  a  flow  drift  tube  (FDT)  at  298  K 
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over  a  wide  range  of  center-of-mass  kinetic  energies  at  low 
pressures.  A  rate  constant  of  6  X 10““  cm^  s“'  has  been  re¬ 
ported  by  Lindinger  at  zero  field  (0.038  eV),  which  decreases 
to  a  minimum  at  0.25  eV  kinetic  energy  and  then  increases  at 
higher  kinetic  energy.  The  rate  constant  has  been  observed  to 
be  independent  of  pressure  over  the  limited  pressure  range 
0.139-0.253  Torr  studied.  Hiraoka®  measured  the  rate  con¬ 
stant  for  reaction  (1)  from  64  to  552  K  using  a  high  pressure 
mass  spectrometer  (HPMS)  at  ca.  3  Torr  and  found  that  the 
rate  constant  decreases  monotonically  with  increasing  tem¬ 
perature  over  the  entire  temperature  range.  Matsuoka  et  al. 
have  also  found  a  total  reaction  rate  constant  at  294  K  for 
reaction  (1)  of7.lXlO““  cm^  s“  *  using  a  time-resolved  at¬ 
mospheric  pressure  ionization  mass  spectrometer  (TRAPI).*® 
There  is  good  agreement  among  the  room  temperature  ther¬ 
mal  rate  constants  reported  by  Smith  et  al,  Lindinger,  and 
Matsuoka  et  al 

Product  formation  has  also  been  studied.  Smith  et  al 
have  determined  that  the  product  ions  formed  in  reaction  (1) 
are  70%  NO'*'  and  30%  NOj  at  300  K.’  By  contrast,  the  drift 
tube  results  of  McCrumb  and  Wameck"  and  the  TRAPI 
measurements  of  Matsuoka  et  al^^  indicate  that  the  majority 
product  is  NO2  .  Dunkin  et  al^  suggest  that  an  additional 
“switching”  reaction  channel  results  in  a  third  product  chan¬ 
nel  giving  O2  ,  proceeding  through  an  intermediate  excited 
state  of  a  (N'*'  -02)*  complex  that  can  decay  into  NO'''  +  0 
and  +  N  products,  or  can  be  stabilized  collisionally  or 
radiatively  to  form  NO^ .  Matsuoka  et  al have  also  re¬ 
ported  that  the  NO2  ions  formed  are  rapidly  lost  in  collisions 
with  N2  at  atmospheric  pressure  with  an  effective  decay  con¬ 
stant  of  3.0Xl0''s“',  carefully  ruling  out  any  interfering 
reactions  of  NO2  .  Hiraoka®  measured  the  relative  branching 
ratios  of  the  NO2  and  NO'*'  products  from  64  to  552  K  and 
found  that  [NO2  ]/[NO''']  decreases  monotonically  with  in¬ 
creasing  temperature  which  is  in  general  agreement  with  the 
values  of  Matsuoka  et  al  fi-om  303  to  384  K.‘^ 

Matsuoka  et  al  further  addressed  the  problem  of  ther¬ 
mal  stability  of  the  NOj  product  of  reaction  (1).  They  deter¬ 
mined  that  the  decay  of  NOj  into  NO'*'  has  a  strong  tem¬ 
perature  dependenee,  increasing  exponentially  with 
increasing  temperature.  Furthermore,  they  have  found  that 
NO2  from  reaction  (1)  is  reactive  and  thermally  dissociates, 
while  the  NO2  formed  in  the  reaction  of  Oj  with  NO2  is 
thermally  stable  and  unreactive.  Therefore,  they  conclude 
that  the  two  NOj  ions  do  not  have  the  same  structure.  These 
authors  further  argue  that  the  atomic  arrangement  in  the 
NO^  from  reaction  (1)  should  be  the  higher  energy  NOO'*' 
linear  isomer.'^ 

However,  the  HPMS  measurements  of  Hiraoka  indicate 
that  ONO"*  is  the  isomer  formed  by  reaction  (1),  based  on 
the  lack  of  thermal  unimolecular  dissociation  of  NO2  at 
465.7  K  and  ab  initio  calculations  of  the  successive  binding 
energies  of  N2  ligands  in  NOj(N2)n  clusters  containing 
ONO'*' .  The  calculated  energies  agree  with  the  measured  val¬ 
ues,  whereas  the  preliminary  calculations  of  Hiraoka  and 
Yamabe  show  that  NOO'*'  is  highly  unstable.’^  Differences  in 
the  nucleophilic  properties  of  N2  vs  O2  further  imply  that  the 
observed  product  is  ONO'*'.® 

In  order  to  elucidate  the  kinetics  of  the  N^  -4-  O2  system. 


Popovic  et  al. 

the  rate  constants  and  product  ion  branching  ratios  for  reac¬ 
tion  (1)  were  measured  fi'om  120  to  1400  K  using  two 
complementary  instruments.  The  discrepancy  concerning  the 
NO2  product  isomer  was  investigated  by  experimental  and 
theoretical  analyses  of  the  thermal  stability,  reactivity,  ener¬ 
getics,  and  structure  of  the  various  isomeric  fonns  of  NO2 
(including  barriers  to  rearrangement).  The  combination  of  a 
SIFT  and  high-temperature  flowing  afterglow  (HTFA)  allows 
this  wide  temperature  range  to  be  probed.  Comparing  the 
current  rate  constants  to  the  previous  measurements  as  a 
function  of  temperature®  and  kinetic  energy*  shows  the  ef¬ 
fects  of  different  types  of  energies  on  the  reactivity.  A  poten¬ 
tial  energy  profile  has  been  calculated  to  elucidate  the  reac¬ 
tion  mechanism. 

II.  EXPERIMENT/COMPUTATION 
A.  Kinetics  measurements 

An  experimental  study  of  the  N^  +  O2  reaction  was  per¬ 
formed  using  two  separate  techniques  in  order  to  encompass 
the  temperature  range  from  120  to  1400  K.  The  SIFT  has 
been  used  in  the  temperature  range  between  120  and  500  K, 
while  the  HTFA  measurements  have  spanned  the  350-1400 
K  range.  Details  of  the  two  systems  are  described 
elsewhere.®**''^  All  of  the  reagent  gases  have  been  used  with¬ 
out  further  purification. 

In  the  SIFT,  N^  ions  are  generated  fi'om  N2  in  an  effu¬ 
sive  source  in  a  high-pressure  source  region.  High  purity 
molecular  nitrogen  is  taken  from  the  boil-off  of  a  high- 
pressure  liquid  nitrogen  dewax.  Two  mechanisms  for  N^ 
generation  are  operative:  three-body  association  of  N'*'  with 
N2  and  Penning  ionization  of  N2 ,  the  latter  of  which  is 
shown  in  Eq.  (2), 

N* +N2-^N3+-tN.  (2) 

Mass  selection  of  N^  using  a  quadrupole  mass  filter  before 
injection  into  the  fast  flow  of  He  buffer  gas  (AGA,  99.995%) 
produces  N^  in  the  flow  tube  in  >99%  purity.  After  reaction, 
the  ions  are  sampled  through  a  blunt  nose  cone  aperture  and 
are  analyzed  by  a  quadrupole  mass  analyzer,  and  then  de¬ 
tected  by  an  electron  multiplier.  Rate  constant  measurements 
are  performed  at  two  reaction  distances,  36  and  52  cm,  and 
the  data  do  not  vary  between  the  two  reaction  lengths  within 
the  experimental  uncertainty.  Measurement  of  the  exponen¬ 
tial  decay  of  the  primary  ion  signal  as  a  function  of  the  O2 
concentration  (Massachusetts  Oxygen,  99.999%)  under 
pseudo-first-order  conditions,  coupled  with  the  previously 
measured  ion  velocities,  yields  the  rate  constant.  The  average 
relative  error  in  the  rate  constants  is  ±15%  with  absolute 
errors  of  ±25%.'^ 

The  product  ion  branching  ratios  are  determined  in  the 
SIFT  by  measuring  the  relative  magnitude  of  the  NOj  and 
NO'*'  signals  obtained  from  extrapolating  to  zero  O2  concen¬ 
tration: 

%[N02®']=  100X[N02"']/([N02+]  +  [NO'*']).  (3) 

In  addition,  the  ratio  of  NO2  to  NO'*',  r,  is  obtained  as 
shown  in  Eq.  (4), 
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r  =  %[NO^]/(100-%[NO2"']).  (4) 

This  ratio  has  been  measured  over  a  relatively  high  range  of 
oxygen  concentrations,  approaching  0.2%  of  the  total  gas 
flow.  No  secondary  chemistry  is  observed  and  the  percentage 
of  NO^  is  independent  of  O2  flow.  The  relative  uncertainty 
of  the  branching  ratios  is  ±10%  of  the  major  product 
peaks.’’ 

In  the  HTFA,  the  ions  are  generated  in  an  ion  source 
region  positioned  orthogonal  to  the  quartz  reaction  tube.*'* 
The  source  has  a  0.5  in.  diameter  exit  aperture  for  these 
experiments,  so  that  the  source  pressure  is  higher  than  that 
found  in  the  flow  tube.  Nitrogen  is  introduced  downstream 
from  a  thoriated  iridium  filament  biased  to  100-150  V  with 
0.1-4  mA  emission  current.  The  same  two  reaction  pro¬ 
cesses  generate  N^  ions  as  in  the  SIFT.  A  helium  buffer  flow 
up  to  30  000  SCCM  is  introduced  through  the  source  region 
and  carries  the  ions  into  the  flow  tube  where  a  commercial 
furnace  controls  the  temperature  (SCCM  denotes  cubic  cen¬ 
timeter  per  minute  at  STP).  Typical  operating  pressures  in 
the  flow  tube  are  between  0.8  and  1.5  Torr.  O2  is  added  52.7 
cm  upstream  from  the  nose  cone  sampling  plate  in  the  region 
where  the  ions  have  reached  thermal  equilibrium.  A  small 
fraction  of  the  flow  is  again  sampled  and  the  ions  are  ana¬ 
lyzed  by  a  quadrupole  mass  analyzer,  and  then  detected  by 
an  electron  multiplier.  The  Nj’^  ion  decay  rate  does  not 
change  with  the  N2  source  gas  flow  at  rates  of  125-500 
SCCM,  indicating  that  the  production  of  the  ions  is  complete 
before  the  reactant  inlet.  Residual  water  contaminant  in  the 
helium  buffer  is  reduced  using  a  liquid  nitrogen  cooled  sieve 
trap. 

At  higher  temperatures,  alkali  ions  themiionically  emit¬ 
ted  from  the  quartz  flow  tube  are  inevitable,  but  these  ions 
are  nonreaetive  background  species  in  the  mass  spectrum. 
Invariance  of  the  alkali  ion  signal  as  a  function  of  O2  flow 
has  been  used  as  an  indicator  that  the  flow  tube  chemistry 
has  properly  equilibrated  as  demonstrated  previously.'*  At 
1400  K,  the  level  of  the  '"K’*'  signal  was  comparable  to  the 
level  of  the  N^  signal  {mlz  =  A2),  and  the  mass  resolution  of 
the  quadrupole  filter  was  adjusted  so  that  the  two  mass  peaks 
are  separated.  Average  relative  errors  of  the  N^  deeay  rate 
are  ±12%  with  absolute  errors  of  ±25%.'® 

Unfortunately,  a  pure  N^  signal  is  not  possible  in  the 
HTFA  because  of  the  lack  of  upstream  mass  selection  and  the 
meehanism  for  N^  ion  generation.  Due  to  the  complexity  of 
the  reactions  involving  N''',  N^ ,  N^ ,  N4  ,  and  NO"’’  also 
present  in  the  flow  tube,  it  is  impossible  to  determine  the 
branching  ratios  for  the  two  reaction  channels  by  simulta¬ 
neously  following  the  two  products,  especially  since  NO’*'  is 
also  a  product  of  a  number  of  background  reactions,  e.g., 
N'''-I-02.  Fortunately,  none  of  the  contaminant  ions  react 
with  O2  to  produce  an  ion  at  mlz  =  A6,  the  mass  of  NO^ . 
Therefore,  the  NO2  fraction  at  a  given  eoncentration  of  oxy¬ 
gen  has  been  evaluated  in  the  HTFA  from  Eq.  (5)  below, 

%[NO+]  =  100([NO+]-[NO+]o)/([N+]o-[N3+]), 

(5) 

where  the  differences  between  the  ion  signals  at  a  given  O2 
concentration  and  zero  O2  eoncentration  are  used.  The  rela- 


FIG.  1.  Temperature  dependence  of  the  rate  constants  for  the  reaction  of  Nj" 
with  O2  as  recently  measured  in  the  selected  ion  flow  tube  (SIFT)  and 
high-temperature  flowing  afterglow  (HTFA).  Also  included  are  measure¬ 
ments  performed  with  a  flow  drift  tube  (Ref  8),  TRAPI  (Ref  10),  SIFT 
(Ref.  7),  flowing  afterglow  (FA)  (Ref  6),  and  HPMS  (Ref  9).  The  Langevin 
collision  rate  constant  for  this  reaction  is  7.0X  10“''  cm^  s“'. 

tive  uncertainty  in  the  product  distributions  determined  this 
way  is  also  10%).  This  uncertainty  is  determined  from  the 
agreement  with  the  SIFT  results  shown  below. 

B.  Ab  initio  calculations 

Energies  of  the  various  isomers  of  the  NO2  speeies  on 
the  ground  state  singlet  potential  energy  surface  have  been 
calculated  at  different  levels  of  theory.  Here  we 
present  only  the  energetics  of  the  NOO'*'  isomer  optimized 
using  multireference  eonfiguration  interaction  at  the 
MRCISD(Q)(15e,12o)/cc-pVTZ  level.  The  detailed  theoreti¬ 
cal  results  will  be  presented  in  a  separate  paper. 

A  potential  energy  profile  has  also  been  calculated.  Ge¬ 
ometries  were  optimized  using  density  funetional  theory 
(DFT)  at  the  Becke3-Lee,  Yang,  and  Parr  (B3LYP)/aug-cc- 
pVTZ  level,  where  aug-cc-pVTZ  denotes  Dunning’s  correla¬ 
tion  consistent  polarized  valence  triple-zeta  basis  set  aug¬ 
mented  with  diffuse  functions."*  Frequency  calculations  have 
been  performed  at  the  same  level  as  the  geometry  optimiza¬ 
tion  to  determine  whether  the  species  are  stable  molecules  or 
transition  states.  Intrinsic  reaction  coordinate  (IRC)  calcula¬ 
tions  have  been  carried  out  to  confirm  the  connectivity  of 
every  optimized  transition  state  at  the  same  level.  To  obtain 
better  energetics,  single  point  energy  calculations  using 
coupled  cluster  theory  at  the  CCSD(T)/aug-cc-pVQZ  level 
have  been  performed  on  the  B3LYP/aug-cc-pVTZ  optimized 
geometries,  where  aug-cc-pVQZ  denotes  Dunning’s  correla¬ 
tion  consistent  polarized  valence  quadmple  zeta  basis  set 
augmented  with  diffuse  functions.'^ 

III.  RESULTS  AND  DISCUSSION 
A.  Rate  constants 

Figure  1  shows  the  rate  constants  of  the  N^  +  O2  reac¬ 
tion  as  a  function  of  temperature  as  measured  in  the  SIFT 
and  HTFA  covering  the  120-1400  K  range;  error  bars  denote 
our  ±15%  experimental  uncertainty.  The  rate  constants  de¬ 
crease  substantially  with  temperature  up  to  1200  K,  reaching 
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FIG.  2.  Dependence  of  the  rate  constants  for  the  reaction  of  Nj  with  O2  on 
the  sum  of  the  average  total  energy  (translational,  rotational,  and  vibrational) 
for  the  present  data  and  those  obtained  in  the  FDT  experiment  (Ref.  8). 


Temperature  (K) 


FIG.  3.  Temperature  dependence  of  the  product  ion  distribution  in  the  reac¬ 
tion  of  N3  with  O2 .  The  TRAPI  results  are  from  the  work  of  Matsuoka 
et  al.  (Ref.  12)  and  the  previous  SIFT  results  are  from  the  work  of  Smith 
el  al.  (Ref.  7).  DT  results  are  from  the  work  of  McCrumb  and  Wameck 
(Ref  11)  and  the  HPMS  from  the  work  of  Hiraoka  (Ref  9). 


a  minimum  at  this  temperature  and  then  increase  up  to  1400 
K.  Data  from  the  two  instmments,  the  SIFT  and  the  HTFA, 
agree  in  the  overlapping  temperature  region  within  the  ex¬ 
perimental  uncertainty. 

Measurements  from  the  other  temperature  dependent 
experiments®’^’^’'”  discussed  in  the  Introduction  are  also 
shown  in  Fig.  1  and  agree  well  with  the  present  work,  if  the 
rate  constants  measured  below  ~180  K  in  the  HPMS  are 
excluded.  The  previous  flow  drift  tube  results  are  discussed 
separately  below.  These  values  are  not  shown  because  clus¬ 
ters  of  the  type  (N2)n  are  formed  and  because  the  rate 
constants  do  not  correspond  purely  to  as  the  reactant  ion. 
Below  180  K,  the  HPMS  rate  constants  show  a  steeper  in¬ 
crease  with  decreasing  temperature,  indicating  that  the 
(N2)n  ions  react  faster  than  does  bare  This  interfer¬ 
ence  is  not  a  problem  in  the  current  experiments  because  no 
N2  is  present  in  the  flow  tube  of  the  SIFT  and  clustering  to 

does  not  occur  at  high  temperature  in  the  HTFA.  The 
room  temperature  value  measured  by  Matsuoka  et  al.  at  at¬ 
mospheric  pressure  is  slightly  larger  than  the  rest  of  the  mea¬ 
surements  and  may  have  been  slightly  influenced  by 
clustering.'”  However,  the  atmospheric  pressure  value  does 
agree  within  the  experimental  error.  All  of  the  rate  constants 
shown  are  considerably  smaller  than  the  Langevin  collision 
rate  constant  of  7.0X  10“ '”  cm^  s“ ' . 

Also  shown  in  Fig.  1  are  rate  constants  from  a  flow  drift 
tube  experiment  (FDT)  taken  as  a  function  of  center-of-mass 
kinetic  energy.*  The  kinetic  energy  is  converted  to  a  transla¬ 
tional  temperature  by  {KE^j^)  =  2il2kT^.  The  agreement 
with  the  drift  tube  data  is  good  at  low  temperature,  but  at 
temperatures  around  800  K  and  higher,  the  data  diverge  for 
reasons  to  be  discussed  below.  In  addition,  the  current  results 
suggest  that  the  rate  constants  begin  to  increase  after  reach¬ 
ing  a  minimum  at  a  translational  temperature  lower  than  that 
found  in  the  drift  tube.* 

In  order  to  compare  the  role  of  various  types  of  energies, 
the  data  are  replotted  as  a  function  of  average  total  energy  in 
Fig.  2.  The  rotational  and  vibrational  temperatures  of  the 
reactant  ions  in  the  FDT  experiments  are  determined  from 
the  center-of-mass  collision  energy  for  ion-buffer  gas  colli¬ 
sions  ((A'£'(,uffer)  =  3/2A:rrot)-^*^  Note  that  a  N2  buffer  was 


used  in  the  drift  tube  experiments.  The  average  vibrational 
energy  is  calculated  as  an  ensemble  average  over  a  Boltz¬ 
mann  distribution  of  harmonic  vibrational  frequencies  of  the 
reactants.®  At  low  temperatures,  vibrations  predomi¬ 
nantly  contribute,  while  the  O2  vibrations  contribute  only  at 
temperatures  over  about  1000  K.® 

The  rate  constants  from  all  of  the  data  sets  decrease  with 
increasing  energy  and  reach  a  small  plateau  in  the  0.2-0.4 
eV  range  before  decreasing  again  to  a  minimum  at  0.7  eV. 
Above  this  energy,  the  rate  constants  increase  with  increas¬ 
ing  energy.  The  temperature  and  kinetic  energy  data  are  in¬ 
distinguishable,  indicating  that  all  forms  of  energy  that  are 
excited  in  this  range  control  reactivity  approximately  the 
same.  Of  course,  cancellation  of  effects  is  possible,  but  past 
experience  indicates  this  occurrence  is  unlikely.^'  The 
equivalency  of  translational  and  rotational  energies  is  a  typi¬ 
cal  behavior  for  most  ion-molecule  reactions.  Vibrational  ex¬ 
citation  often  has  a  different  effect  than  the  other  forms  of 
energy.^' 

B.  Branching  ratios 

Previous  results  on  the  branching  ratios  for  reaction  (1) 
are  few  in  number  and  contradictory.  A  SIFT  experiment  by 
Smith  et  al.  yields  a  value  of  r  =  [N02  ]/[NO‘'']  =  0.43  at 
300  K.^  However,  a  preference  for  reaction  channel  (lb)  (r 
=  1.33)  was  observed  in  the  drift  tube  measurement  by 
McCrumb  and  Wameck  at  the  same  temperature."  Matsuoka 
et  al.  reported  that  r  decreases  from  r=  1.72  at  303  K  to  r 
=  1  at  384  K  in  a  TRAPI  experiment,  which  suggests  that 
reaction  (lb)  is  preferred  at  lower  temperatures  and  that  this 
preference  may  be  reversed  at  higher  temperatures,  i.e.,  r 
<  1.*^  Thus,  the  results  from  an  earlier  SIFT  experiment’  are 
the  opposite  of  the  results  from  the  drift  tube"  and  the 
TRAPI  experiment.'’  The  present  work  attempts  to  (a)  ad¬ 
dress  the  discrepancy  in  the  value  of  r,  (b)  verify  its  tempera¬ 
ture  dependence,  and  (c)  elucidate  the  nature  of  the  NO2 
product. 

Branching  ratios  as  a  function  of  temperature  obtained 
from  the  current  experiments  are  shown  in  Fig.  3,  given  as 
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FIG.  4.  Potential  energy  profile  of  the  NO2  singlet  iso¬ 
mers  optimized  at  the  MRCISD(Q)/cc-pVTZ  level  of 
theory.  On  the  structures,  numbers  are  bond  distances 
and  angles  (in  angstrom  and  degree).  For  the  energetics, 
plain  numbers  are  with  zero-point  energy  (ZPE)  correc¬ 
tion,  and  numbers  in  the  parentheses  are  without  ZPE 
correction.  Only  ONO+_disTS2  was  calculated  at  the 
MRCISD(Q)/cc-pVTZ//CASSCF/cc-pVTZ  level. 


the  ratio  r  defined  in  Eq.  (4).  The  most  striking  feature  of  the 
present  results  is  a  monotonic  exponential  decrease  of  r  with 
increasing  temperature  over  several  orders  of  magnitude,  i.e., 
a  substantial  decrease  in  the  amount  NO2  observed  for  reac¬ 
tion  (1)  at  high  temperature.  This  downward  trend  is  in  ex¬ 
cellent  agreement  with  the  HPMS  temperature  dependent 
data  up  to  ca.  450  K,  indicating  that  the  crossover  to  r<  1 
predicted  by  Matsuoka'^  does  occur.  The  present  data  further 
indicate  that  the  trend  continues  until  the  highest  tempera¬ 
tures.  There  are  no  obvious  explanations  for  (a)  the  exponen¬ 
tial  nature  of  the  decline  and  (b)  for  the  disagreement  of  the 
highest  temperature  HPMS  point. 

The  other  anomalous  point  on  the  branching  ratio  plot  is 
from  the  earlier  SIFT  study.^  In  contrast  to  the  present  study 
where  is  injected  cleanly  into  a  helium  buffer,  the  previ¬ 
ous  SIFT  value  was  obtained  by  injecting  N'*’  into  a  mixed 
He  (0.45  Torr)/N2  (0.2  Torr)  buffer  and  letting  N'*'  cluster 
with  N2  to  form  .  A  potential  problem  would  occur  if 
some  N’''  had  survived  into  the  reaction  zone  (but  not  nec¬ 
essarily  to  the  mass  spectrometer)  since  the  reaction  of  N’’’ 
with  O2  produces  but  not  NO2  In  addition,  Mat- 
suoka  et  a/.'°  have  shown  that  the  NO^  fonned  in  reaction 
(1)  subsequently  reacts  with  N2  with  a  rate  constant  of  1.2 
X  10“  cm^  s“ '  to  form  NO'*' .  This  reaction  would  partially 
convert  NOj  to  NO"''  in  the  previous  SIFT  experiment  and 
skew  the  branching  ratio  to  a  value  less  than  1 . 

C.  NO2  structures  and  unimolecular  reactions 

The  structure  of  the  NO2  ions  produced  in  reaction  (1) 
has  also  been  the  subject  of  debate.  Matsuoka  and 
co-workers'°’*^  have  reported  that  the  NO^  ions  fonned  by 
reaction  (lb)  in  the  TRAPI  are  lost  in  collisions  with  N2  at 
atmospheric  pressure  with  an  effective  decay  constant  of 
3.0X10'*s“'  or  a  bimolecular  rate  constant  of  ~3 
X  10“ cm^  s“ ' .  Furthennore,  they  observed  that  this  decay 
(a)  produces  NO'*',  (b)  is  independent  of  the  O2  and  N2 
content  in  their  reaction  gas  mixture,  and  (c)  is  temperature 
dependent.  In  their  earlier  paper,*®  they  assigned  the  reactiv¬ 
ity  to 

N0j-HN2-^N0+-hN20  (6) 


with  a  bimolecular  rate  of  1.2X  10“'^ cm®  s“'.  In  their  sub¬ 
sequent  paper,  they  attribute  the  loss  of  NO2  to  thermal  de¬ 
composition  of  an  isomer  having  a  NOO®'  structure,  mainly 
because  the  rate  did  not  change  with  buffer  gas  density.  The 
observed  decay  rate  of  NO2  measured  over  a  relatively  nar¬ 
row  range  between  294  and  384  K  using  TRAPI  increases 
exponentially  with  increasing  temperature,  giving  a  rate  of 
(5.6X  10’)exp(-4500/i?r)  s“’.  Furthermore,  they  have 
found  that  NO^  from  reaction  (1)  reacts  strongly  with  C2H6 
to  form  C2H^  .  However,  the  NO2  formed  in  the  reaction  of 
O2  with  NO2  does  not  react  with  C2H6.  Therefore,  they 
coneluded  that  the  two  NO2  ions  do  not  have  the  same  struc¬ 
ture,  also  supporting  the  thermal  decomposition  hypothesis. 

In  contrast,  the  HPMS  measurements  of  Hiraoka  indicate 
that  ONO'*'  is  the  isomer  formed.®  Loss  of  NO2  is  not  ob¬ 
served  at  465.7  K  in  the  HPMS.  Ab  initio  calculations  of  the 
successive  binding  energies  of  N2  ligands  in  NO2  (N2)  clus¬ 
ters  with  ONO'*'  agree  with  the  measured  values,  whereas 
their  preliminary  calculations  show  that  NOO'*'  is  highly 
unstable.’®  In  the  experiments  of  Matsuoka  et  al.  at  atmo¬ 
spheric  pressure,’®’’®  the  single  exponential  decay  of  the 
NO2  signal  indicates  that  only  one  isomer  is  formed. 

To  address  this  issue,  both  high-level  theoretical  calcu¬ 
lations  and  an  experimental  test  of  the  reactivity  of  the  NO^ 
produced  in  reaction  (1)  have  been  performed.  The  details  of 
the  calculations  are  presented  in  a  separate  paper,  with  only 
the  most  relevant  results  at  the  highest  level  of  calculation, 
MRCISD(Q)/cc-pVTZ,  presented  here  in  Fig.  4.  The  calcu¬ 
lations  show  that  the  NOO®'  isomer  is  stable,  with  an  energy 
of  4.51  eV  above  the  ground  state  linear  ONO®'  isomer.  The 
NOO'*'(’2''')  isomer  is  linear  with  a  N-0  bond  distance  of 
1.146  A  and  an  0-0  bond  distance  of  1.233  A.  A  barrier  of 
0.89  eV  prevents  the  isomer  from  dissociating  to 
NO'*'(’S''')  +  0(’Z)),  which  is  only  0.060  eV  endothermic. 
Production  of  NO'’’(’S''')  and  0(®P)  from  NOO‘''(’S'*')  is 
exothermic  by  almost  2  eV,  correlating  to  three  triplet  sur¬ 
faces  accessed  by  curve  crossings.  Two  minima  of  the 
singlet-triplet  seam-of-crossing  (MSX)  have  been  located: 
MSX_’2;'’'/®n  ca.  0.07  eV  below  NOO+_TS  and 
MSX_’S'’'/®X  ca.  0.02  eV  below  NOO''‘_TS,  as  shown  in 
Fig.  4.  Spin-orbit  calculations  indicate  that  both  of  them  have 
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FIG.  5.  Singlet  potential  energy  surface  for  the  reaction 
of  +  O2  calculated  at  the  CCSD(T)/aug-cc-pVQZ// 
B3LYP/aug-cc-pVTZ  level  of  theory.  On  the  structures, 
numbers  are  bond  distances  and  angles  (in  angstrom 
and  degree).  For  the  energetics,  ZPE  correction  has 
been  applied. 


very  strong  interactions:  65.5  cm’"'  for  each  of  the  com¬ 
ponents  with  and  35.7  cm”'  for  the  MSX_'S‘'’/^S, 
suggesting  that  predissociation  could  easily  occur.  On  the 
singlet  ground  state,  there  is  also  a  bent  CAi)  isomer  of 
ONO''"  with  an  ZONO  of  80°.  This  bent  fonn  lies  4.03  eV 
above  the  most  stable  linear  ONO'*’ ,  and  a  small  barrier  of 
0.23  eV  relative  to  the  bent  isomer  links  them.  This  barrier  is 
relatively  small  and  thermal  isomerization  in  the  presence  of 
buffer  gas  should  be  fast  enough  to  prevent  this  species  from 
being  observed  in  our  apparatus.^^  However,  the  barrier  be¬ 
tween  linear  NOO'''  and  bent  ONO'*‘  isomers,  NOO'''_isoTS, 
is  fotmd  to  be  quite  large,  1.36  eV  from  the  linear  NOO'*' 
isomer.  The  dissociation  processes  of  ONO'*^  isomers  are 
also  investigated.  Two  corresponding  transition  states, 
ONO‘^_disTSl  and  ONO'''_disTS2,  connecting  the  linear 
and  bent  isomers  are  found  to  be  4.77  and  1.01  eV  above  the 
corresponding  isomers,  respectively.  Note  that  the  activation 
energy  found  by  Matsuoka  et  a/.'®’'*  for  NO2  reacting  with 
N2  is  much  smaller  than  the  barrier  found  here  for  dissocia¬ 
tion  to  NO"^,  indicating  that  the  process  observed  is  probably 
not  thermal  dissociation.  (Note  also  that  using  the  new  ener¬ 
getics  for  N3,  the  production  of  N2  +  0('D)-I-N0’''  is  0.99 
eV  exothermic.)  Therefore,  all  of  the  transition  states  and 
isomers  except  the  NOO'''_isoTS  discussed  above  lie  below 
the  entrance  channel  energy  for  N^  -f  O2  at  300  K  as  shown 
in  Fig.  4. 

It  is  relatively  easy  to  determine  if  a  high  energy  form  of 
NOj  is  generated  in  the  present  experiments.  A  sufficiently 
large  amount  of  O2  can  be  added  upstream  in  the  flow  tube 
of  the  SIFT  to  drive  reaction  (1)  to  completion  before  a  sec¬ 
ond  gas  is  added  downstream  to  test  the  chemistry  of  the 
product  ions.  What  is  required  is  a  gas  that  does  not  react 
with  the  ground  state  ONO’*'  isomer,  but  will  react  with  the 
NOO"''  isomer  (assmning  that  the  bent  ONO‘'‘  species,  if 
formed,  thermally  isomerizes  to  the  linear  ground  state).  Xe¬ 
non  is  a  good  candidate  since  it  is  a  generally  unreactive  gas. 
For  this  test,  Xe  (Air  Products,  99.995%)  is  introduced  into 
the  SIFT  at  a  distance  30  cm  downstream  from  the  O2  inlet. 


In  a  separate  experiment  for  comparison  purposes,  ground 
state  ONO'*'  is  generated  in  the  source  region  of  the  SIFT  by 
electron  impact  on  nitrogen  dioxide  (Matheson,  99.5%). 
These  ONO"'  ions  are  mass-selected  and  then  reacted  with 
Xe  in  the  normal  manner.  NO2  injected  into  the  flow  tube 
from  the  ionization  of  nitrogen  dioxide  in  the  source  does  not 
react  with  Xe.  By  contrast,  the  46  amu  ion  generated  by 
reaction  (1)  is  totally  consumed  when  Xe  is  added  to  the 
SIFT.  The  reaction  turns  out  to  be  more  complicated  than 
expected,  giving  two  products: 

N00’^  +  Xe^N0'''  +  (Xe+0)  +  356  kJmoF'  (6a) 

->XeO+  +  NO+108  kJmor'.  (6b) 

The  bulk  of  the  reaction  produces  NO^.  However,  the  reac¬ 
tion  is  sufficiently  exothermic  so  that  the  neutral  product 
could  be  either  XeO  or  Xe-fO.  This  difference  in  reactivity 
confirms  that  the  NO2  made  in  reaction  (1)  is  a  high-energy 
isomer  and  is  not  in  the  linear  ONO'*'  form,  in  agreement 
with  the  TRAPI  results.'^  Production  of  XeO’*'  is  rare,  but  it 
has  been  observed  previously.^''  We  have  also  confirmed  that 
C2H6  reacts  with  NOO'''  from  reaction  (1)  to  produce  C2H^ 
and  that  linear  ONO'*'  does  not. 

D.  Potential  energy  surface  and  reaction  mechanism 

Lindinger  suggests  that  the  decrease  of  the  rate  constant 
with  kinetic  energy  (or  temperature)  is  due  to  the  lifetime  of 
the  complex  intermediate  decreasing  with  the  increasing 
relative  velocity  of  the  reactants.*  This  assumes  that  product 
formation  is  slow  compared  to  that  for  dissociation  into  re¬ 
actants  and  that  no  barrier  is  above  the  energy  of  the  reac¬ 
tants.  This  assertion  is  probably  qualitatively  correct  but 
somewhat  simplistic  for  a  quantitative  description.  The  in¬ 
crease  at  higher  temperatures/energies  is  probably  attribut¬ 
able  to  a  new  mechanism  that  becomes  available. 

In  order  to  elucidate  the  probable  mechanisms,  CCSD(T) 
calculations  of  the  potential  energy  profile  have  been  made 
as  mentioned  in  the  Experiment/Computation.  The  resulting 
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FIG.  6.  Triplet  potential  energy  surfaee  for  the  reaction 
of  +  O2  calculated  at  the  CCSD(T)/aug-cc-pVQZ// 
B3LYP/aug-cc-pVTZ  level  of  theory.  On  the  structures, 
numbers  are  bond  distances  and  angles  (in  angstroms 
and  degree).  For  the  energetics,  ZPE  correction  has 
been  applied. 


singlet  and  triplet  potential  surfaces  are  shown  in  Figs.  5  and 
6,  respectively.  The  reactants  are  triplets  and  correlate  to  sin¬ 
glet  and  triplet  products.  The  connectivity  of  the  transition 
states  has  also  been  confinned  as  discussed  in  the 
Experiment/Computation. 

On  the  singlet  surface  shown  in  Fig.  5,  two  intermediates 
are  formed  directly  from  the  reactants.  The  first  intermediate 
(INT_1)  has  the  O2  axis  nearly  parallel  to  the  quasilinear 
(ZN-N-N=168°)  axis,  with  ZN-N-O=105°  and  ZN-0- 
0=115°.  This  “trans”  structure  lies  0.51  eV  below  the  en¬ 
ergy  of  the  reactants.  The  other  “cis”  intermediate  (INT_2)  is 
more  v  shaped  and  lies  0.56  eV  below  the  reactants.  A  tran¬ 
sition  state  (TS_4)  separates  the  two  intermediates  and  is 
1.20  eV  above  INT_1.  INT_1  correlates  only  to  the  NOO''' 
product,  but  the  reaction  must  pass  through  TS_1  with  a  bar¬ 
rier  0.17  eV  above  the  reactants.  Therefore,  it  is  not  impor¬ 
tant  at  low  energies  or  temperatures.  INT_2  is  cormected  to 
both  NO’^  and  NOO'*'  through  separate  transition  states 
(TS_3  and  TS_2)  that  are  0.15  and  0.01  eV,  respectively, 
above  the  energy  of  1NT_2. 

On  the  triplet  surface  shown  in  Fig.  6,  a  V-shaped  inter¬ 
mediate  (INT_3)  is  0.54  eV  below  the  reactants  and  connects 
to  NO'^  and  N2O  singlet  products  through  transition  state 
(TS_6)  only  0.04  eV  above  INT_3.  However,  INT_3  can  only 
be  reached  through  a  cyclic  N30^  transition  state  complex 
(TS_5)  that  is  2.09  eV  higher  in  energy.  Another  transition 
state  (TS_7)  lies  0.53  eV  above  the  reactants  and  simulta¬ 
neously  produces  three  fragments:  N2,  and 

oCp). 

The  features  of  the  surfaces  illustrated  in  Figs.  5  and  6  at 
least  qualitatively  explain  the  reactivity.  At  low  temperatures, 
the  reaction  occurs  at  an  appreciable  fraction  of  the  collision 
rate  and  decreases  with  increasing  energy  and  temperature. 
That  trend  is  typical  for  ion-molecule  reactions  involving 
surfaces  with  stable  intermediates  followed  by  barriers  to 
reaction.  The  simplest  explanation  is  that  decay  of  the  inter¬ 
mediates  into  reactants  involves  a  loose  transition  state  and 


into  products  involves  a  tight  one.  This  picture  is  a  more 
complete  description  of  the  mechanism  than  proposed  by 
Lindinger. 

At  low  energies,  reaction  (1)  occurs  on  the  singlet  sur¬ 
face  in  Fig.  5.  Once  1NT_2  forms,  it  can  decay  back  into 
reactants  or  into  either  product  channel,  NOO'*’  or  NO"^. 
NOO'*’  has  the  lower  barrier  and  is  thus  favored  at  low  en¬ 
ergy.  However,  the  barrier  to  NO'*'  production  is  not  much 
higher  and  the  data  show  that  NO"^  is  favored  by  400  K.  Any 
molecules  that  form  intermediate  INT_1  must  dissociate  back 
to  reactants  unless  there  is  enough  energy  to  overcome  the 
barrier.  However,  sufficient  energy  is  not  available  until 
-800  K  (-0.7  eV). 

The  upturn  in  the  rate  constants  occurs  at  a  total  energy 
of  —0.7  eV,  Fig.  2.  Two  reaction  paths  may  contribute  to  this 
trend  at  higher  energies.  First,  TS_4  on  the  single  surface 
connecting  the  two  intennediates  INT_1  and  INT_2  has  a 
barrier  of  0.69  eV,  approximately  coinciding  with  the  energy 
at  the  upturn.  Therefore,  this  additional  mechanism  leading 
to  products  may  begin  to  contribute  to  the  overall  reactivity 
above  0.7  eV,  causing  an  increase  in  the  total  rate  constant. 
Second,  reaction  also  becomes  possible  on  the  triplet  surface 
in  Fig.  6  where  the  0.54  eV  baixier  to  TS_7  can  be  overcome, 
opening  up  the  channel  that  produces  three  fragments  includ¬ 
ing  NO'*'.  In  addition,  a  further  reduction  in  the  NOO”^  prod¬ 
uct  at  higher  energies  can  come  from  unimolecular  decom¬ 
position  of  NOC)'''  giving  NO'^.  The  multireference 
configuration  interaction  calculation  results  shown  in  Fig.  4 
indicate  that  a  crossing  from  the  singlet  NOO’*'  surface  to  a 
triplet  surface  at  MSX_*S'''/^n  occurs  about  0.8  eV  above 
NOO'''  minimum  that  allows  predissociation  to  NO''^('S''') 
+  0(^P)  that  is  2  eV  lower  in  energy.  This  pathway  is  con¬ 
sistent  with  the  NOO'*'  unimolecular  decomposition  results 
of  Guilhaus  et  al.  showing  that  NOO'''  dissociation  produces 
NO'*'  fragments  with  2  eV  kinetic  energy.^^  The  present  cal¬ 
culations  also  indicate  that  a  barrier  of  0.89  eV  exists  to 
dissociate  NOO"*  into  NO'''  and  0(*D).  Since  the  products 


9438  J.  Chem.  Phys.,  Vol.  121,  No.  19,  15  November  2004 

can  be  formed  hot  (1.1  eV  exothermic),  decomposition  can 
happen  at  relatively  lower  temperatures,  in  addition  to  the 
mechanisms  just  discussed. 

IV.  CONCLUSIONS 

The  kinetics  of  the  reaction  of  with  O2  has  been 
studied  over  a  wide  temperature  range,  up  to  1400  K.  The 
rate  constants  are  found  to  decrease  from  120  to  1200  K 
before  increasing  up  to  1400  K  and  agree  with  previous  mea¬ 
surements  taken  over  limited  temperature  ranges.  Comparing 
the  rate  constants  to  those  obtained  in  a  drift  tube  shows  that 
all  forms  of  energy  behave  similarly. 

Two  products  are  observed,  namely,  NO2  and  NO'*' .  The 
former  is  favored  at  low  temperature  and  the  latter  at  higher 
temperature.  The  ratio  r  of  NOj  to  NO^  decreases  exponen¬ 
tially  with  temperature  from  r  =  2  at  120  K  to  r<0.01  at 
1400  K.  In  this  paper,  a  straightforward  test  of  the  reactivity 
on  the  NO2  product  ion  was  made  that  shows  conclusively 
that  this  product  is  a  high  energy  isomer,  most  likely  NOO"^ . 
Theoretical  calculations  confirm  that  this  is  a  stable  isomer, 
with  a  barrier  to  dissociation  of  about  0.9  eV.  The  energy  of 
the  NOO'*'  isomer  is  determined  to  be  4.5  eV  above  that  of 
the  most  stable  isomer,  linear  ONO’*' .  The  calculated  poten¬ 
tial  energy  surfaces  for  the  reaction  giving  singlet  products 
explain  most  of  the  current  data. 
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